clc;
clear;
f=@(x,y)exp(-x.^2/2).*sin(x.^2+y);
tic,
[x,y,z]=intfunc2(f,-2,2,-1,1);%逐步积分结果存入z
toc,
surf(x,y,z),
I=z(end,end)%z中取最后一个元素作为积分结果
